Calculating the trap density of states in organic field-effect transistors from 
experiment: A comparison of different methods 
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The spectral density of localized states in the band gap of pentacene (trap DOS) was determined 
with a pentacene-based thin-film transistor from measurements of the temperature dependence 
and gate-voltage dependence of the contact-corrected field-effect conductivity. Several analytical 
methods to calculate the trap DOS from the measured data were used to clarify, if the different 
methods lead to comparable results. We also used computer simulations to further test the results 
from the analytical methods. Most methods predict a trap DOS close to the valence band edge 
that can be very well approximated by a single exponential function with a slope in the range of 
50 — 60meV and a trap density at the valence band edge of ~ 2 x 10 21 eV -1 cm -3 . Interestingly, 
the trap DOS is always slightly steeper than exponential. An important finding is that the choice 
of the method to calculate the trap DOS from the measured data can have a considerable effect on 
the final result. We identify two specific simplifying assumptions that lead to significant errors in 
the trap DOS. The temperature-dependence of the band mobility should generally not be neglected. 
Moreover, the assumption of a constant effective accumulation layer thickness leads to a significant 
underestimation of the slope of the trap DOS. 

PACS numbers: 73.61. Ph, 73.20.At, 73.20.Hb, 85.30.De 

Keywords: organic semiconductor, thin-film transistor, trap density of states, transistor modeling 
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I. INTRODUCTION 

World-wide research on organic field-effect transistors 
is at a high level as this new technology is poised to 
enter the market. [J0 An appealing feature of this tech- 
nology is that organic semiconductors can be deposited 
by thermal evaporation or from solution at low cost on 
large areas while keeping the substrates close to room 
temperature. Consequently, organic semiconductors are 
promising candidates for future flexible and low-cost elec- 
tronics. 

The mobility of charge carriers in organic field-effect 
transistors is comparable to the mobility in hydrogenated 
amorphous silicon thin-film transistors (lcm 2 /Vs) and 
thus is already adequate for many applications. [3, U IE [|| 
In addition to a high mobility, useful organic transistors 
must have a near zero threshold voltage, a steep sub- 
threshold swing and a high electrical and environmental 
stability. The transistor parameters and stability of or- 
ganic field-effect transistors are intimately related to the 
efficiency of the charge transport mechanism and the ex- 
tend of charge carrier trapping in extrinsic traps. The 
main scientific challenge thus is to clarify the nature of 
the charge transport mechanism and the microscopic ori- 
gin of charge carrier traps in organic field-effect transis- 
tors. 

Field-effect transistors can be used to determine the 
underlying spectral density of localized states in the 
band gap, i.e. the trap densities as a function of 
energy (trap DOS). This has been extensively done 



'Electronic address: kalb@phys.ethz.ch 



with thin-film transistors (TFT's) employing amorphous 
semiconductors or with TFT's based on polycrystalline 

silicon. 0, 1 M M El m m M HE M Ea M, 51 This 

approach is expected to be of great value for the under- 
standing of any novel semiconductor in a field-effect tran- 
sistor including organic small molecule semiconductors, 
polymeric semiconductors or ZnO.[2fl| 

The research efforts to calculate the trap DOS from 
measurements of organic field-effect transistors have in- 
creased only recently. On the one hand, the density of 
states function can be calculated from the linear regime 
transfer characteristics in a st raig htforward fashion with 
an analytical methodjp, H HE HI [H, [H S3, H HI 
H [M 1H [H, [H, Hi] This approach has the advan- 
tage of giving an unambiguous result but errors may re- 
sult from the various simplifying assumptions. On the 
other hand, a density of states function can be postu- 
lated a priori and the corresponding transistor charac- 
teristics can be calculated by means of a suitable com- 
puter program. The density of states function is then 
iteratively refined until good agreement between the mea- 
sured characteristic and the computer-simulated curve is 
achieved. [M HI [H, [M IS EL El 

Several analytical methods have been used to calcu- 
late the trap DOS from transfer characteristics of organic 
field-effect transistors. In order to eventually clarify the 
microscopic origin of charge carrier traps in organic field- 
effect transistors, it is highly desirable to quantitatively 
compare the results from various experiments and from 
various research groups. An important prerequisite to 
such a systematic comparison is to test if the different 
analytical methods lead to comparable results. We ap- 
plied the most frequently used analytical methods to the 
same set of measured transfer characteristics. Moreover, 
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we applied a computer simulation program to determine 
a trap DOS that leads to simulated transfer characteris- 
tics closely matching the measured data. 

In the following we are dealing with pentacene and p- 
type conduction. For convenience the charge carriers are 
called holes although they may be of pronounced pola- 
ronic character. [43j Moreover, we use terms such as va- 
lence band edge, band mobility or effective density of ex- 
tended states. However in order to apply the calculation 
methods described below, we do not necessarily need to 
have band transport or the existence of extended states. 
The calculation methods may be applied as long as the 
charge transport can be described by a transport level 
with a distribution of localized states below this trans- 
port level (trap-controlled transport).}^ 

We begin by summarizing the widely-used analytical 
description of an organic field-effect transistor. This de- 
scription is only valid for samples with a low trap density 
and negligible contact resistances. The trapping and re- 
lease times are assumed to be much shorter than the time 
necessary to measure a transistor characteristic, i.e. we 
have no current hysteresis. In Sec. IIIII we present the 
trap DOS calculated with the different methods for a 
pentacene TFT with a Si02 gate dielectric, i.e. a sample 
with a significantly high trap density. The essential equa- 
tions of the different analytical methods are also given in 
this section along with specific details about the use of 
the calculation methods. 



The drain current in the saturation regime (\Vd\ > 
\V g — Vt\) quadratically depends on gate voltage, i.e. 



Id = 



{V g -V t f 



(2) 



Fitting a straight line to the square root of the mea- 
sured drain current yields the band mobility (1q and the 
threshold voltage Vt- This ideal behaviour can be ob- 
served in organic ficld-cffcct transistors with a low trap 
density. For example in Fig. [T] we show the near-ideal 
transfer characteristic of a pentacene single crystal ficld- 
cffcct transistor (SC-FET) with a Cytop™ fiuoropoly- 
mer gate dielectric. (28j 

The onset voltage V on and the subthreshold swing S 
are other important device parameters. The onset volt- 
age is defined as the gate voltage where the drain current 
exceeds the noise level which typically is at 10 -12 A (see 
Fig. [5]). The subthreshold swing is a measure of how eas- 
ily a transistor can be switched from the off-state to the 
on-state. It is denned as [45[ 



S = 



dV a 



d(logJd)' 



(3) 



With the simplistic assumption that both the density of 
deep bulk traps Ni> u ik and the density of (deep) interface 
traps Ni n t are independent of energy, the subthreshold 
swing may be written as [4(j| 



fcTlnlO 



II. ANALYTICAL DESCRIPTION OF AN 
IDEAL FIELD-EFFECT TRANSISTOR 

In the linear regime (|Vd| <C \V g — Vt\), the drain current 
Id of an ideal field-effect transistor is given by 



W 

id = Y^Ci(v g - v t )v d , 



(i) 



where Cj is the capacitance per unit area, W and L are 
the channel width and length, V g and Vd are the gate 
voltage and the drain voltage and is the band mobility 
which is independent of gate voltage. Eq.[T] predicts a lin- 
ear dependence of the drain current on the effective gate 
voltage Id oc (V g — Vt). A linear regression of the mea- 
sured transfer characteristic thus yields the band mobility 
Ho and the threshold voltage Vt- The threshold voltage 
is defined as the gate voltage above which essentially all 
of the incrementally added gate-induced charge is mo- 
bile ( "free" ) . The threshold voltage depends on the trap 
density in the device and on the flatband voltage Vfb- 
The flatband voltage is the gate voltage which needs to 
be applied in order to enforce flat bands at the insulator- 
semiconductor interface. A non-zero flatband voltage can 
result from a difference of the Fermi level in the semicon- 
ductor and in the gate electrode. More importantly, the 
flatband voltage is influenced by charge that is perma- 
nently trapped at the insulator-semiconductor interface 
or within the gate dielectric. 



1 + — ( V tsN bu ik + eN int ) 



This may be simplified as follows: 



S 



fcTlnlO 



(4) 



(5) 



Both the deep bulk traps and the interface traps con- 
tribute to the trap density Njj (per unit area and unit 
energy). [U HgJ The subthreshold swing thus is a sim- 
ple measure of the deep trap density. Fig. [2] shows the 
same data as in Fig. Q] on a logarithmic scale. The sub- 
threshold swing is as steep as S = 0.3V/dec. With 
Eq. [5] and Ci — 4.3nF/cm 2 we calculate a trap density 
from S as low as 2Vp = 1.1 x 10 11 cm~ 2 eV _1 . Assuming 
an effective accumulation layer thickness of a = 7.5 nm 
this results in a volume trap density of N = N^/a = 
1.5 x 10 17 cmr 3 eV -:L . 

If the experimental transfer characteristics are linear in 
the linear regime and quadratic in the saturation regime, 
the approach described above is self-consistent and the 
extracted mobility /io and threshold voltage Vt have a 
clear meaning, i.e. (j,q is the band mobility and Vt is volt- 
age above which the incrementally added gate-induced 
charge is placed in the valence band. In samples with 
an increased trap density, the drain current in the linear 
regime may however increase faster than linearly. For 
example, this can be seen in the lower panel of Fig. [3] 
where we show the measured transfer characteristics of 
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FIG. 1: (Color online) High-performance pentacene single- 
crystal field-effect transistor (SC-FET) with a Cytop™ flu- 
oropolymer gate dielectric. The graph shows the square root 
of the drain current (full black line) and the dashed red line is 
a linear fit of the measured data. The square root of the drain 
current linearly depends on gate voltage in accordance with 
the well-known field-effect transistor equation for the satu- 
ration regime. This linear dependence is a mark of the low 
trap density at the insulator-semiconductor interface, as well 
as the high field-effect mobility of fi = 1.4cm 2 /Vs and the 
near- zero threshold voltage. 
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FIG. 2: (Color online) Same data as in Fig. [TJ of a pen- 
tacene/Cytop SC-FET plotted on a logarithmic scale. The 
forward and reverse sweeps are shown. The very low trap 
density in the active region of the transistor manifests itself 
in several desirable properties: near-zero onset voltage V on , 
very steep subthreshold swing of 5* = 0.3 V/dec, as estimated 
from the dashed red line, and negligible current hysteresis. 



a pentacene TFT with SiC>2 gate dielectric for several 
temperatures. The transconductance {did/ dV g )v d now 
increases monotonically with gate voltage. The percent- 
age of the gate-induced holes that arc free increases with 
gate voltage and this leads to the "superlinear" trans- 
fer characteristics. Strictly speaking, the threshold volt- 
age is not reached even at relatively large gate volt- 
ages. Eq. [1] and Eq. [5] are not suitable for this type of 
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FIG. 3: (Color online) Transfer characteristics of a pentacene- 
based thin-film transistor (TFT) with a Si02 gate dielec- 
tric measured at various temperatures (symbols). The up- 
per and lower panel show the same data on a logarithmic 
and linear scale. The red lines are computer-simulated trans- 
fer characteristics. The increased subthreshold swing (upper 
panel) and the reduced field-effect mobility (jj, ~ 0.2cm 2 /Vs 
at V g = —50 V and T = 299 K) are a result of a relatively high 
trap density. The high trap density also results in the drain 
current to increase faster than linearly in the linear regime 
(lower panel). 

transistor. [49l |50| This is similar to the case of amorphous 
silicon field-effect transistors, where the trap densities are 
substantial. [5l| We note that for the TFT, the gate ca- 
pacitance is Ci = 13.3nF/cm 2 . This is about three times 
larger than the capacitance for the SC-FET in Fig. Q] and 
Fig. [2] The difference between the two types of transis- 
tors is even more drastic than it appears when comparing 
the graphs. 

III. QUANTIFICATION OF THE TRAP DOS 

Field-effect transistors can be used to quantify the un- 
derlying trap DOS if a more sophisticated description of 
the device physics is used. We applied several analytical 
methods to the linear-regime transfer characteristics in 
Fig. [3] We also applied a simulation program to deter- 
mine a trap DOS that leads to simulated transfer charac- 
teristics closely matching these measured transfer char- 
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acteristics. Moreover, the results were compared to the 
crude estimate of the trap density from the subthreshold 
swing (Eq. We used a dielectric constant of e.; = 3.9 
for SiC>2 and e s — 3.0 for pentacene.[12, HH The thick- 
ness of the SiC>2 gate dielectric was I = 260 nm and the 
pcntaccne film was d — 50 nm thick. The channel was 
L — 450 fim long and W — 1000 /im wide. 

All analytical methods and the simulation program are 
based on the following simplifying assumptions: 

- the organic semiconductor is homogenous perpendicu- 
lar to the insulator-semiconductor interface, and 

- insulator surface states only introduce an initial band 
bending without applied field, i.e. contribute to a non- 
zero fiatband voltage Vfb- 

As a consequence we obtain an effective trap DOS. In 
the case of TFT's with polycrystalline films, the trap 
densities to be determined are an average over intra-grain 
and inter-grain regions and may also be influenced, to 
some extend, by trap states on the surface of the gate 
dielectric. 



A. Analytical methods 

Several additional assumptions are made to simplify 
the analytical methods: 

- the charge density is homogenous along the transistor 
channel (from source to drain), and 

- for the trapped holes, the Fermi function is approxi- 
mated by a step function (zero temperature approxi- 
mation), and 

- the valence band is approximated as a discrete energy 
level at the valence band edge Ey with an effective den- 
sity of extended states Ny. The occupation of these 
extended states is calculated with the Boltzmann func- 
tion, and 

- the temperature dependence of the Fermi energy Ep as 
well as of the interface potential Vo is neglected (neglect 
of the statistical shift). 

The first assumption is appropriate only if the trans- 
fer characteristics are measured at a low drain voltage. 
In that case we can assume the "unperturbed" situa- 
tion where charge is accumulated by a gate voltage in 
a metal-insulator-semiconductor (MIS) structure but no 
drain voltage is applied. 50] 

The final results from the different methods are sum- 
marized in Fig. 2] and Table [J along with the result from 
the simulation program. The parameters Nq and Eq in 
Table U were obtained by fitting to each trap DOS an 
exponential function 

N(E) = N exp(-E/E ). (6) 



TABLE I: Parameters resulting from different methods to cal- 
culate the trap DOS and the trap DOS from each method is 
shown in Fig. [4] The different methods were applied to the 
same set of measured data (Fig. [3} • To obtain the parameters 
iVo and Eo, an exponential function N(E) = No exp(—E/Eo) 
was fitted to the calculated trap DOS in each case. Most 
methods lead to a slope in the range of Eo = 50 — 60meV 
and to a trap density at the valence band edge of iVo ~ 
2 x 10 21 eV _1 cm -3 . Some methods also lead to an estimate 
of the band mobility po- 



Method 


Ref. 


No 

(eV _1 cm- 3 ) 


Eo 

(meV) 


Mo" 

(cm 2 /Vs) 


Kalb I 


[32] 


8.5 x 10 21 


60 


0.7 


Kalb II 


This paper 


2.4 x 10 21 


59 




Fortunato 


[18, 26] 


1.1 x 10 21 


60 




Griinewald 


[8,IQ] 


4.0 x 10 21 


41 




Computer sim. 


[40] 


1.5 x 10 21 


50 


0.3 


Lang 6 


[25] 


1.1 x 10 20 


115 


17 


Horowitz 


M 


7.0 x 10 20 


48 


0.2 


Subthreshold 
swing 6 c 


[46, 47,48] 


N = 4.4 x 10 is 





"Band mobility at T = 299 K. 

^Effective accumulation layer thickness of a = 7.5 nm. 
'Assumption: trap DOS independent of energy. 



Fig. [4] and Table U also contain the trap density as es- 
timated from the subthreshold swing of S — 2.4V/dec 
with Eq. O and an effective accumulation layer thickness 
of a = 7.5 nm. Some methods also lead to an estimate of 
the band mobility /iq (see Table QJ. 

The final results in Fig. 0] and Table U are discussed 
in Sec. IIVI The simulations are described in more detail 
in Sec. IIII B[ In the following we describe the different 
analytical methods and give some specific details about 
how the methods were used. 

All methods require linear regime transfer character- 
istics measured at several temperatures (as in Fig. [3]). 
The exception is the method by Griinewald et al. which 
only requires a single linear regime transfer characteristic 
measured at one temperature (e.g. room temperature). 
In order to calculate the trap DOS from the transfer char- 
acteristics Id{Vg) we need the field-effect conductivity a 
and the field-effect mobility /i at first. Provided that con- 
tact effects are negligible, the drain current in the linear 
regime may be written as 

W 

Id = Y aV d ( ? ) 

and the field-effect conductivity a is 

a = nCi{V g - V FB ). (8) 

Vfb is the fiatband voltage and fi is the (gate-voltage 
dependent) field-effect mobility, i.e. an effective mobil- 
ity in contrast to the band mobility //o- The field-effect 
conductivity can be calculated from 
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FIG. 4: (Color online) Spectral density of localized states in the band gap of pentacene (trap DOS) as calculated with several 
methods from the same set of measured data. The energy is relative to the valence band edge (VB). The estimate from the 
subthreshold swing (dashed grey line) assumes the trap DOS not to depend on energy and can only be regarded as a rough 
estimate for the traps slightly above the Fermi energy. All other methods result in a trap DOS that increases slightly faster 
than exponentially with energy. The choice of the method to calculate the trap DOS has a considerable effect on the final 
result. 



The calculation of the field-effect mobility from measure- 
ments of organic field-effect transistors is still controver- 
sial for samples with an increased trap density. [50] If we 
differentiate Eq. [7] with respect to gate voltage, we have 

M 



did W ( cV 
— = -r^Vd U + (V g - V FB ) — 



dV a 



dV a 



(10) 



'9 ^ \ vv g, 

The field-effect mobility is most often calculated from 



did 

WV d C t V dV, 



(11) 



which means that the second term in Eq. [10] is generally 
neglected. 

Contact effects can introduce significant errors when 
the field-effect conductivity and the field-effect mobility 
are calculated. [HR I5r3| The contact-corrected field-effect 
conductivity can be calculated from gated four-terminal 
measurements according to 



WV' d 



(12) 



L' is the distance between the voltage sensing electrodes 
and V' d = V\ — V 2 is the voltage drop between these 
electrodes. [2S| [32| The effective field-effect mobility fi is 



not influenced by parasitic contact resistances when cal- 
culated from gated four-terminal measurements accord- 
ing to 



V 



WV'Q 



did 
dV n 



(13) 



The mobilities as calculated with Eq. [TTJ or Eq. [Jj2 
overestimate the true field-effect mobilities to some 
extend. [54| This is because for trap-controlled transport, 
the mobility increases with gate voltage. This leads to 
a positive second term in Eq. [10] which is neglected in 
Eq. QT] (and Eq. [BJ. However, the use of Eq. [TT] (or 
Eq. I13[) is advantageous because the definition and ex- 
traction of a flatband voltage is circumvented. More- 
over, this approach to the field-effect mobility is most 
often used which guarantees a good comparability of the 
mobility values. Consequently, we choose Eq. [13] (the 
contact-corrected version of Eq. [TTj) to calculate the field- 
effect mobility in the present study. Alternatively, the 
field-effect mobility could also be calculated from [54[ 



V 



wv>a (y g - v FB ) 



(14) 



In that case we would not need to differentiate the mea- 
sured data but a reliable estimate of the flatband voltage 
(threshold voltage) would be required. [fJij 
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In Fig. [3] we show the measured transfer characteristics 
Id(V g ) from gated four-terminal measurements. In addi- 
tion to the drain current Id(V g ), the potentials Vi (V g ) and 
^2(^9) between the grounded source electrode and the re- 
spective voltage sensing electrode were measured simulta- 
neously while keeping the source-drain voltage constant. 
This was done by connecting the source of the transistor 
to the ground connector of an HP 4155A parameter an- 
alyzer and by measuring the channel potentials V\ and 
V 2 with two additional SMU's in the "source current - 
measure voltage" mode with a sourced current of OA. 
For all analytical methods we used the four-terminal con- 
ductivity a = cr(V g ) as derived from gated four-terminal 
measurements with Eq. [12] and the field-effect mobility 
was calculated according to Eq. [131 This allowed for 
the calculation of a trap DOS that is free from con- 
tact artifacts. [1^, [32j Moreover, we only used currents 
above \IA = 1 nA for the calculation of trap DOS (1 nA 
limit). (HI 

For the following description of the analytical methods, 
Fig. 2 and Fig. 3 in Ref. [32| are useful. 
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FIG. 5: (Color online) Activation energies E a , E' a and E'l 
as determined with linear regressions according to Eq. 1 1 5 1 
Eq. [43] and Eq. 1301 The graph also shows a smooth fit of 
the activation energies in each case (red lines). There is a 
significant difference between E a , E' a and E'^. 



1. Method by Lang et al. 

For this method (Ref. I25I ), the activation energy 
E a (V g ) is defined by 



v(V g )=Ae X pl- 1 £ 



E„ 



kT 



(15) 



and A is assumed to be a constant. The activation en- 
ergy is determined from the measured data with a linear 
regression of In a vs. l/T for each gate voltage according 
to In a = In A — E a /kT. The energetic difference between 
the Fermi level Ef of the sample and the valence band 
edge at the insulator-semiconductor interface is approx- 
imated with the measured activation energy E a (V g ) of 
the field-effect conductivity c, i.e. 



E a ~ Ey — Ep — eV 



(16) 



Vq = \V(x = 0)| is the potential right at the insulator- 
semiconductor interface. The ^-direction is normal to 
this interface. Ey is the energy of the valence band edge 
far from the insulator-semiconductor interface (at x = d, 
Fig. 3 in Ref. I32T ). The underlying idea is the follow- 
ing: a change of the gate voltage by AV g leads to a shift 
of the activation energy E a (i.e. of the effective Fermi 
level Ef — Ef + eVb at the insulator-semiconductor in- 
terface) by AEp ~ AE a . The change in gate voltage 
A.V g corresponds to a total hole density per unit area 
of AP = CiAV g /e. Then, the abrupt approximation is 
made: the charge in the accumulation layer is constant 
up to a distance a from the insulator-semiconductor in- 
terface and zero for larger distances. For the present 
method, it is assumed that the parameter a does not de- 
pend on gate voltage. [25J Consequently, we have a change 
of the volume hole density of Ap — AP/a close to the 



insulator-semiconductor interface. By neglecting the free 
charge we can estimate the trap density to be 



N = A P /AE a 



ea 



AEa 
AV n 



(17) 



If one replaces the difference quotient in Eq. [17] by the 
respective derivative we have the final result 



N(E) 



d ( dE a 



(18) 



Consequently, the function N(E) is calculated with 
Eq. [JS] and N(E) is plotted as a function of the energy 
E = E a (V g ) ^E V -E F - eV M 



The band mobility fio may be estimated with 

'E n 



fi = A* exp 



kT 



(19) 



by introducing the value of the measured activation en- 
ergy E a and the field-effect mobility fi at a fixed and 
sufficiently high gate voltage V^.[57j 

We now give some specific details about the applica- 
tion of this method to our data. The activation energy 
E a = E a (V g ) was determined according to Eq . [TBI and is 
shown in Fig. [5] The activation energy E a (V g ) was then 
represented by a smooth fit (red line in Fig. [5]) in order 
to suppress the noise in the data. We used an effective 
accumulation layer thickness of a = 7.5 nm.(;25] The band 
mobility fj,o was calculated with Eq. [Ill] for a high gate 
voltage of V g = -50 V and T=299K. 



2. Method by Horowitz et al. 



Also for this method (Ref. 21J) the abrupt approxima- 
tion is made. However, in contrast to the method by 
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Lang et al., the present method allows for a gate- voltage 
dependence of the effective accumulation layer thickness 
a = a(V g ). As a consequence of the abrupt approxima- 
tion, the potential V(x) in the organic semiconductor is 
given by 



vw-w(i-i)- 

with the interface potential 
V 



CtU„ 



2e e s 



(20) 



(21) 



U a 



\V g — Vfb\ is the gate voltage above the flatband 



voltage Vfb- With the total hole density per unit area 
P = CiUg and Eq.[2T]we obtain an equation for the total 
volume hole density p which is 



P 

a 



CfU 2 



2e e s eV 



(22) 



Assuming again the abrupt approximation, it can be 
shown that: 



eV = Ey — Ep — fcTln 



(23) 



/i = n(V g ) is the gate-voltage dependent field-effect mo- 
bility as calculated with Eq. [T21 Ny is the effective den- 
sity of extended states. For each temperature, the trap 
DOS is now calculated separately. To do so, a value of 
the product poNy is assumed a priori and the interface 
potential Vq is calculated with Eq. [221 We note that this 
also requires an estimate of the difference between the 
Fermi energy and the energy of the valence band edge 
far from the insulator-semiconductor interface Ey, i.e. 
an estimate of Ey — Ep. Vq = Vo(V g ) from Eq. [221 is 
used to calculate the volume hole density p with Eq. [221 
The Fermi function is approximated by a step function 
(zero-temperature approximation). Its derivative then is 
a delta function. Consequently, the trap DOS is eventu- 
ally obtained by numerically differentiating the hole den- 
sity p from Eq. [22] with respect to the interface potential 
V from Eq. [221 i.e. 



e aVo 



(24) 



The trap densities from Eq. [2H are finally plotted as a 
function of the energy E — Ey — Ep — eVb as calculated 
with Eq. [531 The trap DOS from the measurements at 
different temperatures will generally not coincide at first. 
The procedure is thus repeated for different values of the 
product noNy until the trap DOS curves calculated from 
the data taken at different temperatures, coincide. 

The band mobility po is calculated from the final pa- 
rameter noNy by assuming a value of the effective density 
of extended states Ny. If Ny is fixed, po is the adjustable 
parameter. 
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FIG. 6: (Color online) Trap densities calculated with the 
method by Horowitz et al. from the measurements at dif- 
ferent temperatures. The band mobility p,o is an adjustable 
parameter in this method. Upper panel: po = 2cm 2 /Vs, 
lower panel: po = 0.2 cm 2 /Vs. For po = 0.2cm 2 /Vs the trap 
densities from the measurements at different temperatures co- 
incide. 



We proceed by giving specific details about the use of 
this method. We estimated that Ey — Ep = 0.5 eV. In 
order to determine U g — \V g — Vfb\, the flatband volt- 
age Vfb was taken to be equal to the device onset volt- 
age at room temperature. We thus have Vfb = —4.6 V. 
Moreover, the effective density of extended states Ny 
was assumed to be equal to the density of the pentacene 
molecules, i.e. Ny = 3 x 10 21 cm~ 3 . For example, the 
trap densities were calculated from the measurements at 
different temperatures with /iq = 2cm 2 /Vs and the re- 
sult is shown in the upper panel of Fig. The trap 
densities from the different temperatures do not coin- 
cide. The procedure was repeated for different values of 
Ho- For (iq — 0.2cm 2 /Vs the energetic distributions of 
traps do coincide (lower panel of Fig. [5]). This trap DOS 
was taken as the final result along with the band mobility 
of Mo = 0.2cm 2 /Vs. 



Method by Fortunato et al. 



This method (also called temperature method) is de- 
scribed in Ref. Ila and [H The trap DOS is calculated 
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with 



N(E) 



if 



2e dV a 2 



dV(x 



dx 



(25) 



The electric field dV/dx in Eq 
dV(x) 



| may be written as 



dx 



e t U g - V 



x=0 



I 



(26) 



The interface potential Vq in Eq. [23 and Eq. is ob- 
tained as described in the following. First of all, the 
derivative of the field-effect conductivity can be written 
as 



da 
dVa 



Mo 



N v e e l 
lp(Vo) 



exp 



Ey — Ep — eV 
kT 



(27) 



and both the band mobility /x and the exponential factor 
depend on temperature. [18| It can be shown, that the to- 
tal hole density p{Vq) varies much less with temperature 
than the exponential term in Eq. [27j [HI As described in 
the following, the calculation of a normalized field-effect 
conductivity a' eliminates the temperature dependence 
due to the band mobility fM in Eq. [27] [H| The field- 
effect mobility at RT and at each reduced temperature T 
are calculated with Eq. [T3J Then, the normalized field- 
effect conductivity a' is calculated for each temperature 
according to 



RT 



T ' 

Mo 



(28) 



,RT 



/i'" and /j, in Eq. [28] are the field-effect mobilities at 
RT and at a reduced temperature T evaluated at a fixed 
and sufficiently high gate voltage. /j, rt and /j,q are the 
respective band mobilities. [l|J[26[ Clearly, the derivative 
of the normalized field-effect conductivity may now be 
written as 



da' 
dVa 



/V r7777\ ex P 



lp(Vo) 



Ev — Ef 



eV 



kT 



(29) 



and the exponential term is the only term with a temper- 
ature dependence (the temperature dependence of p(Vo) 
is neglected). Consequently, the activation energy E" as 
determined with linear regressions according to 



da' ( El 

cx exp - — - 



dVa 



kT 



(30) 



is approximately equal to the difference between the 
Fermi energy and the valence band edge at the insulator- 
semiconductor interface, i.e. 



E" ^E V -E F ~ eV . 



(31) 



Once E"(V g ) is known, the interface potential Vo(V g ) can 
be calculated with Eq. [3T] assuming a priori a value for 
the difference Ey — Ep. Then we can calculate the elec- 
tric field with Eq. [221 The result is finally introduced 



in Eq. [25] and the numerical differentiation with respect 
to Vo from Eq. [3T] eventually yields the trap DOS as a 
function of energy E = E" w E v — Ep — eV . 

It is instructive to consider the statistical shift in this 
context. Both the Fermi energy Ep and the interface 
potential Vo depend on temperature. For the moment 
we assume that this temperature-dependence is linear, 



Ep — Ev, 



Ep + aT 



and 



eV = eV ° + (a - (3)T. 



(32) 



(33) 



a and (3 are constants and E F and V® are the Fermi 
energy and the interface potential at T — OK. From 
Eq. [SH and Eq. [33] we see that Ey — Ep — eV = 
Ey - Ep- eV$ + 0T. If this is introduced into Eq. [29] 
we have, within this linear approximation, a constant 
prefactor exp(/3//c) and the activation energy i?" is in 
fact a better approximation of the difference between the 
Fermi energy and the valence band edge at the insulator- 
semiconductor interface at T = OK, i.e. 



K ~E V -E%- eV °. 



(34) 



Nevertheless, the energy scale is not corrected for the 
statistical shift for the present method. The trap DOS is 
simply plotted as a function of E". Neglecting the sta- 
tistical shift is a common feature of all analytical meth- 
ods in the present comparison that employ temperature- 
dependent measurements (all methods except for the 
method by Griinewald et al.). 

Again, we assumed that Ey — Ep = 0.5 eV and Vfb = 
—4.6 V. In order to determine the activation energy 
with Eq. [3D] the field-effect mobilities were calculated 
with Eq. [T3J leading to a mobility that increases mono- 
tonically with gate voltage at all temperatures. These 
functions were evaluated at V g — —50 V. The mobili- 
ties at 299, 289, 280, 270, 259 and 249 K respectively 
are fi = 0.17, 0.16, 0.15, 0.13, 0.11 and 0.09 cm 2 /Vs. ^ 
E" is significantly different from the activation energy 
E a (Fig. . Again, the activation energy was repre- 
sented by a smooth fit in order to suppress the noise in 
the data (red line in Fig. [5]). The smoothed function was 
used for the calculation of the trap DOS. 



4- Method by Griinewald et al. 

This method does not require temperature-dependent 
measurements. H [13, 0) HH R allows to convert a sin- 
gle transfer characteristic into the underlying density of 
states function which may be advantageous in certain sit- 
uations where temperature-dependent measurements are 
not possible. [5^| In addition, it is not necessary to con- 
sider the temperature dependence of the Fermi energy, 
the interface potential or the band mobility. Moreover, 
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the method is not based on the abrupt approximation. 
The interface potential Vo as a function of gate voltage 
is calculated from 



exp 



kT e s lao 



eVo 
kT 

u, 



U g a(U g ) 



' kT 
a(U g )dU g 



1 



(i 



(35) 



For each gate voltage, Eq. [35] is numerically evaluated 
using the measured field-effect conductivity a (Eq. [T2"]) . 
Eventually, we have the complete function Vq = Vo(V g ). 
The total hole density p can be calculated with Vo ac- 
cording to 



P(V ) 



J 2 e 9 V dU n 



dV Q 



(36) 



Within the zero-temperature approximation the trap 
DOS N(E) can then be written as 



N(E) 



1 dp{Vo) 
e dVo ' 



(37) 



This means that we do a numerical differentiation of the 
hole density from Eq. [36] with respect to the interface 
potential Vq from Eq. 0J5] The trap DOS can be plotted 
as a function of the energy E = eVb, i.e. as a function of 
the energy relative to the Fermi energy Ep. It can also 
be plotted as a function of the energy E = Ey — Ep — eVb 
which requires the difference Ey — Ep to be estimated. 

For the present method, the gated four-terminal mea- 
surement at T — 299 K was considered. Again, we as- 
sumed E v -Ep = 0.5 eV and V FB = -4.6V. 



Method I by Kalb et al. 



The free hole density P/ ree per unit area is written as 
P free « a(V g )N v exp ( E v-^-eV^ (gg) 



with 



a(V g 



kT 



2kTe Q e s 



2m -1 eQU n 



(39) 



a(V g ) is the effective thickness of the accumulation layer. 
E = kT is the slope of the trap DOS and m = T /T.\j$ 
Since the field-effect conductivity can be written as 



(T = efioPfre 



(40) 



the difference Ey — Ep — eVb is approximated by the ac- 
tivation energy E a (V g ) of the field-effect conductivity a. 
E a is determined with linear regressions from the mea- 
sured data according to Eq. 1151 This procedure implies, 
that the temperature dependence of the mobility (j,q as 



well as the temperature dependence of the effective ac- 
cumulation layer thickness a are negligible compared to 
the exponential temperature dependence. By substitut- 
ing dVo — —dE a /e in Eq. 1371 and Eq. [M] we finally have 
the trap DOS 



" IK, 



dEa 

dU a 



(41) 



as afunction of the energy E — E a (V g ) ~ Ey — Ep — eVb- 
The band mobility fiQ can be estimated with 

fi = cr/{eP free ), (42) 

where a is the measured field-effect conductivity and 
Pf re e is calculated according to Eq. 



32] 



We give some specific details about the use of this 
method: We used V FB = -4.6 V. The trap DOS was 
calculated with Eq. [IT] from the smooth fit of the acti- 
vation energy E a in Fig. [5] The parameter m = Tq/T 
in Eq. [3S] is only relevant for the calculation of the band 
mobility [1q with Eq. 1421 To obtain this parameter, an 
exponential function N(E) = Noex.p(—E/Eo) was fit- 
ted to the trap DOS that had previously been obtained 
with Eq. gT] This gave E = kT = 60 meV and thus 
m = 2.33 at T — 299 K. The band mobility was calcu- 
lated for a gate voltage of V g = -50 V and T = 299 K. 



6. Method II by Kalb et al. 

As suggested by Fortunato et al, the temperature de- 
pendence of the band mobility /zo can be eliminated by 
calculating a normalized field-effect conductivity a' for 
each temperature according to Eq. [55] 

In order to improve upon the method by Kalb et al., 
the normalized activation energy E' a (V g ) is determined 
for each gate voltage with a linear regression according 
to 



(43) 



a' in Eq. [43] is the normalized field-effect conductivity 
according to Eq. [55] E' a is a better approximation of 
the difference between the Fermi energy and the valence 
band edge, i.e. we now have 



E„ w E\ 



E f 



eV . 



(44) 



E' a is now used instead of E a in Eq.[4TJ i.e. the trap DOS 
is finally calculated with 



N(E) « — 



eJ 2 



U 



fdE' 



\dU g 



(45) 



It is plotted as a function of the energy E = E' a w Ey 
E F - eVb- 
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7. Influence of the choice of parameters 



We also investigated how the choice the effective ac- 
cumulation layer thickness a and the difference between 
the Fermi level and the valence band edge far from the 
insulator-semiconductor interface Ev—Ef effect the final 
result. The effective accumulation layer thickness a needs 
to be fixed for the method by Lang et al. Clearly, this 
choice significantly affects the trap DOS (Fig. [7]). For the 
method by Horowitz et al, Ey — Ep — 0.5 eV was cho- 
sen. We repeated the calculations for Ey — Ep = 0.8 eV 
and again, the trap densities from the measurements at 
different temperatures were found to coincide with a pa- 
rameter of /in = 0.2 cm 2 /Vs. The results are compared 
in Fig. [Jj The slope of the trap DOS calculated for 
Ey — Ef = 0.5 eV and 0.8 eV are almost identical in 
both cases, but the trap densities are reduced to some 
extend due to the larger value of eVb in the denomi- 
nator of Eq. [22] Also for the method by Fortunato et 
al, the parameter Ey — Ep needs to be known. How- 
ever, this method is not sensitive to the choice of this 
parameter: calculations for Ey — Ep = 0.8 eV give es- 
sentially the same result. Clearly, for the method by 
Giinewald et al., the guess of Ey — Ep has a significant 
influence on the trap DOS since the energy scale is given 
by Ey — Ep — eVo and only eVg is known. Therefore, the 
choice of Ey — Ep = 0.8 eV instead of 0.5 eV leads to a 
parallel shift of the trap DOS by 0.3 eV along the energy 
scale as shown in Fig. [7J For the method by Kalh et al. 
we neither need to make an assumptions about Ey — Ep 
nor about the effective accumulation layer thickness a. 



The necessity to choose the effective density of ex- 
tended states (Ny — 3 x 10 21 cm -3 in the present 
study) leads to an uncertainty in the absolute vale of the 
band mobility. The volume density of extended states 
is given by the density of the pentacene molecules (« 
3 x 10 21 cm" 3 ) multiplied by two due to the spin degree of 
freedom. l59l. [60T| However, a value of as 6 x 10 21 cm -3 
is likely to overestimate the effective density of extended 
states Ny . The underlying spectral density of extended 
states (in cm~ 3 eV -1 ) should drop close to the valence 
band edge. The extended states close to the valence 
band edge are the most important contribution to Ny 
though. Therefore, the volume density of molecules with- 
out the degeneracy factor (Ref. [321 ) or even half the 
molecular density (Ref. [2lh have been used as approx- 
imations of Ny. In the case of pentacene, a value as 
low as Ny = 1 x 10 21 cm" 3 has also been used441] From 
the method by Horowitz et al. it appears, that Ny can- 
not be higher than 3 x 10 21 cm~ 3 . A band mobility of 
= 0.2cm 2 /Vs was calculated from the product ^Ny 
with Ny = 3 x 10 21 cm -3 : a larger Ny would lead to a 
band mobility lower than the field-effect mobility which 
would not be a reasonable result. 
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FIG. 7: (Color online) The choice of certain parameters in 
the analytical methods can have a significant effect on the 
trap DOS. For the method by Lang et al., an effective ac- 
cumulation layer thickness of a = 1.5 nm (dashed red line) 
instead of 7.5 nm (full red line) leaves the slope of the trap 
DOS unchanged but leads to a significant increase of the 
overall trap densities. For the method by Horowitz et al., 
the use of Ev — Ef = 0.8 eV (dashed violet line) instead of 
Ev — Ef = 0.5 eV (full violet line) results in a rather small 
change in the magnitude of the trap densities. The method 
by Grunewald et al. very sensitively depends on the choice of 
Ev — Ef ■ using Ev — Ef = 0.8 eV (dashed green line) instead 
of Ev — Ef =0.5 eV (full green line) results in a parallel shift 
of the trap DOS by 0.3 eV. 



B. Computer simulation of the transfer 
characteristics 

For the modeling of the transfer characteristics, we 
used the Matlab®-based program developed by Ober- 
hoff et a/.[40( The essence of this approach is that the 
program calculates the transfer characteristics from the 
trap DOS, the spectral density of extended states and the 
band mobility fi$. It is assumed that the valence band 
has a rectangular shape, i.e. the density of extended 
states (in cm~ 3 eV~ 1 ) is constant anywhere in the valence 
band. [401 ] The program can simulate the transfer charac- 
teristic at any temperature T as long as the band mobil- 
ity no at this temperature is also fixed a priori. The full 
Fermi-Dirac statistics is included. fioT] This means that, 
in contrast to the analytical methods, the Fermi function 
is not approximated and the temperature-dependence of 
the Fermi energy Ep is not neglected. In Fig. [3] we 
show simulated transfer characteristics (red lines) closely 
matching the measured data (symbols). The trap DOS 
from which these transfer characteristics were calculated 
is also shown in Fig. 0] (light blue line). Table [I] lists the 
respective parameters including the band mobility at RT. 

We now give some more specific details about the sim- 
ulations. The program allows for a consideration of par- 
asitic resistances at the source and drain contacts. For 
the present simulations, we have however assumed negli- 
gible contact resistances. This is supported by the gated 
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four-terminal measurements, which show that the total 
contact resistance is significantly lower than the chan- 
nel resistance at all temperatures. More specifically at 
V g = —50V and T = 299 K, the channel resistance 
is about ten times larger than the contact resistance 
and still five times larger than the contact resistance at 
V g = -50 V and T = 249K. We use 10 22 cm^eV" 1 for 
the density of extended states in the rectangular band. 
In essence, this is an approximation of the density of ex- 
tended states close to the valence band edge since only 
these states are of importance for the charge transport. 
The value is lower by a factor of two compared to the vol- 
ume density of extended states (6 x 10 21 cm -3 ) divided 
by the bandwidth (0.3 eV, Ref. [6l|). The reduced value 
accounts for a drop of the spectral density of extended 
states close to the valence band edge in analogy to the 
choice of Ny — 3 x 10 21 cm" 3 for the analytical methods. 
To obtain a good fit of the transfer characteristics at all 
temperatures, it was necessary to allow for a tempera- 
ture dependence of the band mobility ^o- For the fit in 
Fig. U the band mobilities ^ at 299, 289, 280, 270, 259 
and 249 K were fixed at 0.32, 0.28, 0.25, 0.21, 0.17 and 
0.14 cm 2 / Vs, respectively. 



IV. DISCUSSION 

We begin with a discussion of the results in Fig. 2] and 
Table HI The estimate from the subthreshold swing as- 
sumes the trap densities not to depend on energy. It is 
derived from the subthreshold swing and can thus be re- 
garded as a rough estimate of the density of traps slightly 
above the Fermi energy Ep. It is however gratifying to 
note the agreement between the trap DOS derived from 
the subthreshold swing (valid for E ~ Ep) and the low- 
est values of the trap DOS from the other methods. All 
other methods result in a trap DOS that increases some- 
what faster than exponentially with energy. From Fig. 0] 
and Table U we see, that the choice of the method to cal- 
culate the trap DOS has a considerable effect on the final 
result. 

The method by Kalb et al. gives a good estimate of 
the slope of the trap DOS, but leads to an overestima- 
tion of the overall magnitude of the trap densities, i.e. 
the parameter N . This is because the temperature de- 
pendence of the band mobility fj,o in Eq. [3D] is neglected 
compared to the exponential temperature dependence 
(Eq. |3"8"| . The improved method by Kalb et al. (method 
II) uses the activation energy E' a instead of E a , and E' a 
is calculated from the normalized field-effect conductiv- 
ity at each temperature (Eq. I43[) . This means that the 
temperature dependence of the band mobility is properly 
taken account of. The correction has a considerable ef- 
fect on the overall magnitude of the trap densities. The 
improved method leads to a trap DOS that is in much 
better agreement with the other methods and in partic- 
ular with the result from the simulations. This means 
that the temperature dependence of the band mobility 



Ho should generally not be neglected when calculating 
the trap DOS with an analytical method. Method II by 
Kalb et al. is similar to the method by Fortunato et al. It 
is easier to be used but is based on additional simplifica- 
tions. The difference between the method by Fortunato 
et al. and the method II by Kalb et al. is the use of the 
activation energy E^ of the derivative da' /dV g instead of 
the activation energy E' a of a' . Moreover, Vb in Eq. EBlis 
not neglected by Fortunato et al, contrary to the method 
by Kalb et al. Neglecting Vo does not lead to significant 
differences for sufficiently high gate voltages since the 
interface potential is typically less than 0.5 V. On the 
other hand, from Fig. [5] we see that there are significant 
differences between E' a and E'^. The difference between 
the two methods, therefore, is almost exclusively due to 
the use of E" instead of E' a . The method II by Kalb 
et al. corrects for the temperature-dependence of the 
band mobility but neglects the temperature-dependence 
of a(V g ) (Eq. I39p against the exponential temperature- 
dependence in Eq. [35] This still is a source of error. For 
the present example we have mkT / (2m— 1) = 16.41 meV 
at T = 299,K and mkT/(2m - 1) = 13.07meV at 
T = 249 K, i.e. a ratio of 1.26. This should be com- 
pared to fi RT /n T = 1.9 (fi RT = 0.17 cm 2 /Vs and fi T = 
0.09 cm 2 /Vs at T = 249 K). The method by Fortunato 
et al. is in excellent agreement with the result from the 
computer simulations. 

The method by Lang et al. leads to a significant er- 
ror in the slope of the trap DOS. This is mainly due to 
the assumption of a gate-voltage independent effective 
accumulation layer thickness a. If we allow for a gate- 
voltage dependent effective accumulation layer thickness 
in the context of the abrupt approximation, the effec- 
tive thickness a(V g ) decreases with increasing gate volt- 
age (Eq. |3"9"]) . As a consequence, the assumption of a 
constant thickness a in the denominator of Eq. [18] leads 
to an overestimation of the trap density at low gate volt- 
ages (at energies far from the valence band edge) and 
to an underestimation of the trap density at high gate 
voltages (at energies close to the valence band edge). 

From Fig. [7] we see that some analytical methods do 
not lead to an unambiguous result. In principle, the dif- 
ference Ey — Ep may be approximated with the activa- 
tion energy near the flatband condition. However, this 
activation energy can often not be measured because the 
off-current of an organic field-effect transistor is often due 
to experimental limitations and is not related to the con- 
ductivity of the organic semiconductor. If Ey — Ep is not 
experimentally accessible, we have an uncertainty in the 
trap densities using the methods by Horowitz et al. and 
Griinewald et al. as shown in Fig. [7] Moreover, Fig. [7] 
shows the significant dependence of the result from the 
method by Lang et al. on the choice of the constant ef- 
fective accumulation layer thickness a. The methods by 
Fortunato et al. and by Kalb et al. do not lead to these 
uncertainties. 

The analytical methods approximate the Fermi func- 
tion to a step function for the trapped holes and use 
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Boltzmann's approximation for the free holes. The tem- 
perature dependencies of the Fermi energy Ep and the 
interface potential Vq are also neglected. These assump- 
tions appear to be less restrictive because the trap distri- 
butions from most analytical methods are in good agree- 
ment with the result from the simulations which do not 
involve these assumptions. 

Most analytical methods lead to a band mobility /Lto 
that is comparable to the value of //o from the simu- 
lations. It is important to note that the band mobil- 
ity /Xo from most methods is only slightly higher than 
the field-effect mobility fi at high gate voltages, i.e. 
fi w 0.2cm 2 /Vs at V g = -50 V and T = 299 K for this 
sample. Since \i = {P free I 'P)fJo this means that even in 
samples with an increased trap density (TFT's), most of 
the gate-induced holes are free. 

For the method by Fortunato et al. and for method II 
by Kalb et al., the activation energy (E'£ or E' a ) is cal- 
culated from the normalized field-effect conductivity at 
each temperature. The field-effect conductivity is nor- 
malized to the field-effect mobilities at high gate voltages 
/i RT I n T . This ratio is an approximation of the ratio of 
the respective band mobilities /i RT / (1q . For example, the 
field effect mobility at V g = -50 V is [i RT = 0.17cm 2 /Vs 
at 299 K and [i T = 0.09 cm 2 /Vs at 249 K. This gives 
a ratio of fi RT / /i T = 1.9. [64[ From the simulations we 
have a band mobility of ^ RT = 0.32 cm 2 /Vs at 299 K 
and $ = 0.14 cm 2 /Vs at 249 K. The ratio of the band 
mobilities thus is fi RT / fifi — 2.3. These two ratios are 
very similar, indeed. This further supports the correction 
of the field-effect conductivity as suggested by Fortunato 
et al. 

We also note that from the simulations we have a band 
mobility fio that decreases as the sample is cooled down. 
This may indicate that the trap-free transport process is 
a hopping transport but may also be limited by a ther- 
mally activated process at the grain boundaries. 

Finally, we recall that the trap DOS was calculated 
from transistors with a rather small gate capacitance and 
we thus have relatively large operating voltages. If tran- 
sistors with a high gate capacitance are to be used in 
order to quantify the trap DOS great caution is required. 
Transistors with a sufficiently high gate capacitance can 
be operated at gate voltages of only a few volts (e.g. 
2 — 3 V) . [(52| This is comparable to the magnitude of the 
interface potential Vq (< 0.5 V). However, all methods 
apart from the method by Fortunato et al. assume that 
the total charge per unit area can be approximated ac- 
cording to d ( V g — Vfb — Vq ) w Ci(V g — Vfb)- Significant 
errors are to be expected if Vo is neglected for low- voltage 
operating transistors. The method by Fortunato et al. 
does not neglect Vq and could thus be used in this sce- 



nario. For a low-voltage operating transistor, the choice 
of the parameter Ey — Ep for the method by Fortunato 
et al. is however expected to be a source of ambiguity 
because Vq can no longer be neglected in the numerator 
of Eq. OH 

V. SUMMARY AND CONCLUSIONS 

Several different methods were used to quantify the 
spectral density of localized states in a pentacene-based 
organic thin-film transistor. The trap DOS derived from 
the simple formula for the subthreshold swing is in rather 
good agreement with the lowest values of the trap DOS 
from the other, more sophisticated methods. Most meth- 
ods result in an almost exponential trap DOS close to the 
valence band edge with a typical slope of 50 meV. We find 
that the choice of the method to calculate the trap DOS 
has a considerable effect on the final result. More specif- 
ically, two assumptions lead to significant errors in the 
trap DOS. First, neglecting the temperature dependence 
of the band mobility can lead to a rather large overesti- 
mation of the trap densities. Secondly, the assumption of 
a gate-voltage independent effective accumulation layer 
thickness results in a significant underestimation of the 
slope of the trap DOS. A general conclusion of this study 
is that it is necessary to consider the specific deviations 
of a given calculation method if one compares energetic 
distributions of trap states from organic field-effect tran- 
sistors evaluated by different groups with different meth- 
ods. 

The methods by Fortunato et al. and the method II by 
Kalb et al. do not lead to ambiguities due to the choice 
of parameters, and this constitute a significant advan- 
tage. The computer simulations do not approximate the 
Fermi function and may therefore be seen as the most 
reliable result. Simulating the transfer characteristics at 
various temperatures can, however, be a time consuming 
endeavor due to the large number of possibilities to fix 
the trap DOS and the band mobilities. While all meth- 
ods have their advantages and disadvantages, the method 
by Fortunato et al. is relatively easy to use and gives an 
unambiguous result in excellent agreement with the com- 
puter simulations. 
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